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Anisotropic pairing between fermion species with different fermi momenta opens two-dimensional 
areas of gapless excitations, thus producing a spatiaUy homogeneous state with coexisting superfluid 
and normal fluids. This breached pairing state is stable and robust for arbitrarily small mismatch 
and weak p-wave coupling. 



Introduction. Recently there has been considerable in- 
terest in the possibility of new forms of superfluidity that 
could arise when one has attractive interactions between 
species with different fermi surfaces. This is stimulated 
by experimental developments in cold atom systems 1] 
and by considerations in high-density QCD Possible 
coexistence of superfluidity with gapless excitations is an 
especially important qualitative issue. Spatially homoge- 
neous superfluid states that coexist with gapless modes 
at isolated points or lines in momentum space arise in 
a straightforward way when BCS theory is generalized 
to higher partial waves Gapless states also are well 
known to occur in the presence of magnetic impurities 
and, theoretically, in states with spontaneous break- 
ing of translation symmetry , where the gapless states 
span a two-dimensional fermi surface. Strong coupling 
between different bands also may lead to zeros in quasi- 
particle excitations and gapless states 0. For spheri- 
cally symmetric (s-wave) interactions a spherically sym- 
metric ansatz of this type naturally suggests itself when 
one attempts to pair fermions of two different species 
with distinct fermi surfaces, and a pairing solution can 
be found i, S ES [H. The stability of the result- 
ing state against phase separation or appearance of 
a tachyon in the gauge field (negative squared Meissner 
mass) 13 is delicate, however 0, 0| . It appears to re- 
quire some combination of unequal masses, momentum- 
dependent pairing interactions, and long-range neutral- 
ity constraints. Here we demonstrate another possibility: 
direction-dependent interactions, specifically p-wave in- 
teractions. In this case, stability appears to be quite ro- 
bust. It seems quite reasonable, intuitively, that expand- 
ing an existing (lower-dimensional) locus of zeros into a 
two-dimensional zone should be significantly easier than 
producing a sphere of gapless excitations "from scratch" . 
We shall show that it even occurs for arbitrarily small 
coupling and small Fermi surface mismatch. 

Interactions relevant to pairing can be dominated by p- 
wave (or higher) harmonics under several circumstances. 
If the s-wave interaction is repulsive, it will not be sub- 
ject to the Cooper instability, and will not induce pairing. 
The Cooper instability can be regarded as an enhance- 
ment of the effective interaction for attractive channels 



as one integrates out high-energy modes near the Fermi 
surface. Thus the effective Hamiltonian will come to re- 
semble the form we assume if the interspecies interactions 
are repulsive in s-wave but attractive in p-wave. Fermi 
statistics forbids diagonal (intraspecies) s-wave interac- 
tions; if the higher partial waves are repulsive, or weakly 
attractive, the model discussed here will apply. In the 
context of cold atom systems, tuning to an appropriate 
p-wave Feshbach resonance, as recently reported in [l^ . 
can insure interspecies p-wave dominance. 

A crucial difference between the model we consider 
and the conventional p-wave superfluid system, '^He, lies 
in the distinguishability of the paired species. So al- 
though there are two components, there is no approxi- 
mate quasispin symmetry, and no analogue of the fully 
gapped B phase [l7|. In the absence of a magnetic fleld 
^He has an approximate SO{S)l x S0{3)s x [/(I) sym- 
metry under separate spatial, rotations, spin rotations, 
and number, which is spontaneously broken to the di- 
agonal S0{3)l+s in the B phase. The residual sym- 
metry enforces a gap of uniform magnitude in all di- 
rections in momentum space. The systems we consider 
have quite different symmetry and breaking patterns, for 
example U{1)l, x U{1)a x U{1)b U{1)l,+a+b for 
two spin-polarized species A, B in a magnetic field, or 
S0i3)L X U{1)a X C/(1)b ^ Uil)L^+A+B if the mag- 
netic field can be neglected. The reduced residual sym- 
metry allows for interesting direction-dependent struc- 
ture in momentum space. (In the A phases '^He pairs 
effectively as two separate single-species systems, which 
again is quite different from our set-up.) 

Experimental realizations of p-wave interaction in cold 
atom systems have been reported recently in Ref. |l8j. 
Feshbach resonance in p-wave occurs between '^'^K atoms 
in / = 9/2, nif = —7/2 hyperfine states. This is in con- 
trast to the s-wave resonance, which occurs between non- 
identicdf = 9/2, m/ = -9/2 and / = 9/2, m/ = -7/2 
states 1191. A promising system for possible observa- 
tion of the p-wave breached pairing superconductivity 
is a mixture of / = 9/2, rrif = —9/2 and / = 9/2, 
nif ~ —7/2 atoms ^''K tuned into the repulsive side 
of the s-wave Feshbach resonance. Different densities 
(Fermi momenta) of m/ = —9/2 and mf — —7/2 par- 
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tides can be prepared using different magnitudes of an 
initial additional magnetic field, which is then removed. 
Large atomic relaxation times ensure that the created 
(metastable) states will exist long enough to allow for- 
mation of a superfluid phase. 

Model: Having in mind cold atoms in a magnetic trap 
with atomic spins fully polarized by a magnetic field, we 
consider a model system with the two species of fermions 
having the same Fermi velocity vp, but different Fermi 
momenta pp ± I /vp- The effective Hamiltonian is 

H = ^[6>^ap+e^6lp6_p-A;at,6t_p_Ap6_pap] (1) 
p 
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Ap = ^p-k(ak&_k)- Here the attractive inter-species 
interaction is — Vp-k within the "Debye" energy 2ujo 
around the Fermi surface [up) 3> /), and the intra- 
species interaction is assumed to be either repulsive or 
negligibly small. Excitations of the Hamiltonian 

= ip + A^ + /, are gapless provided that there 

are areas on the Fermi surface where / > Ap. The gap 
equation at zero temperature. 
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can be simplified by taking the integral over d^p, 
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where v — p^^ / {2n'^vp) is the density of states. In the last 
expression it is assumed that / and A„ are dimension- 
less, and scaled with the "Debye" energy: / — > 2ujpl, 
An 2w_DAn. In deriving Eq. © we neglect de- 
pendence of lp_k on the absolute values of p and k; 
this is valid for ujd <C Ep. At weak coupling we 
may linearize in the partial wave expansion, V{n, n') = 
J2i,Tn ^iYhn{n)Yi*^{n'). Assuming p-wave dominance, we 
parameterize V{n, n') = g{n-n') with 5 > 0. P-wave gap 
parameters can arise in the forms Yio and Yi±i, describ- 
ing polar and planar phases respectively. 

Polar phase: An ~ Yio(n). We look for a solution in 
the form An = Acos(n, z) where z is a fixed but arbi- 
trary direction (rotational symmetry is broken) . The gap 
equation becomes. 
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where 9* = arccos^, for z = I/A < 1, and 6* = for 
z > 1. Performing the i nteg rations (detailed calculations 
will be given elsewhere |2flj) we obtain the algebraic gap 
equation. 
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where y = A/Aq is the relative magnitude of the gap 
compared to its value at / = 0, 



Ag°' = exp - 



1.40 exp - 



t^g 
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for weak coupling. There is a factor 3 in the exponent 
with anisotropic interaction instead of 1 as in the s-wave 
BCS. For small values of z the solution to the gap equa- 
tion is y = 1 — 2^ with x — //Aq. We depict the solu- 
tion of the polar phase gap equation in Fig. 1, with the 
following numerical values of the characteristic points: 
XA = (4/37re)i/3 = 0.538, yA = e'^l^ = 0.717 (at the 
point A y'{x) ^ 00), yc = e^^/^ /2 = 0.358 (at the point 
C y{x) = 0). 




FIG. 1: Solutions y{x) of the gap equation in the polar, m = 
0, and planar, m = 1, phases. The lower branch corresponds 
to the unstable state. The branches merge at the points where 
y'(x) — > 00, beyond which there are no non-zero solutions of 
the gap equation. The broken line A = / is included to guide 
the eye. 
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Planar phase: An ~ Yii(n), An ^ Yi_i(n). We now 
look for a solution in the form An — A sin (n, z)e"^, 
where (j) is the polar angle in the plane perpendicular 
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to z. The gap equation becomes, 
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where 6* — arcsinz, for z < 1, and = 7r/2, for z > 1. 
Performing the integration we obtain the algebraic gap 
equation, 



ln(l/y) = -^ + |(3 + z2)ln 



1 + z 
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where again y is the relative magnitude of the gap com- 
pared to its value for a zero mismatch. For the planar 
phase, 
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For small values of z the solution to the gap equation has 
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the form y = 1 ~ and a; is defined as before. Note 
that the planar phase gap is more robust than the polar 
phase, being perturbed by the fourth power instead of the 
third. Solution of the gap equation for the planar phase 
is depicted in Fig. 1, with the following numerical values 
of the characteristic points: xa = 0.674, yA — 0.787, 
ZA = XA/yA = 0.856; xc = e^^/^ = 0.435. 

Stability: The condensation energy is given by (at T — 

0), 
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Evaluating this expression for z = //A < 1, we obtain 
for polar phase 
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which is negative for z < 0.537, and for planar phase 
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which is negative for z < 0.623. For our specific model 
Hamiltonian, at weak coupling, the planar phase is more 
stable. For / > A the condensation energy is always 
positive, indicating that the lower branches are unstable. 

Following the standard methods in the theory of super- 
conductivity |2l| we calculate the super-currents in our 
system under the influence of homogeneous in space vec- 
tor potential A. The super-current is anisotropic, ji — 
^j^KikAk with the components given by {kxx = Hyy), 
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For the polar phase, assuming z > 1 we find. 



1 - 3z3/47r, 

1 - 37rz/4 + 37rz3/ 
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The coefficient k^x becomes negative at z > 0.480 {kzz 
at higher values of z > 0.752) indicating an instability 
with respect to a transition into some inhomogencous 
state (probably similar to a LOFF state). For the planar 
phase, 
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The coefficients Kxx and Kzz remain always positive for 
the whole range of z < 0.623 where the gap equation © 
has stable solutions. Thus, we find that the planar phase 
has lower energy and higher density of Cooper pairs than 
the polar phase and is therefore more stable. 

Specific Heat: The important manifestation of the BCS 
states with gapless excitations is the appearance of the 
term linear in temperature in the specific heat, which 
is characteristic for a normal Fermi liquid. The specific 
heat is given by 
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dT 



dT 



(16) 



where E^ 



+ Al + I. At low temperatures T < 

/ the first term in Eq. (|16ll gives an exponentially small 
contribution. The second term, with E^ , in Eq. H16|l is, 
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Performing the integration, we calculate the contribution 
of the gapless modes to the specific heat at T ^ / to be 
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As expected, the "normal" contribution to the specific 
heat, is proportional to the area occupied by the gapless 
modes, i.e. the I/A strip around the equator for the 
polar phase and the /^/A^ islands around the poles for 
the planar phase. 

Conclusion and Comments: We have presented sub- 
stantial evidence that our simple model supports the 
planar phase gapless superfluidity in the ground state. 
For I <^ A the gapless modes contribute high powers in 
terms of mismatch, ~ for the solution and ~ for 
the heat capacity, i.e. they represent small perturbations. 
The residual continuous symmetry of this state, and its 
favorable energy relative to plausible competitors (nor- 
mal state, polar phase) suggest that it is a true ground 
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state in this model. The planar phase is symmetric un- 
der simultaneous axial rotation and gauge (i.e., phase) 
transformation. Also, we obtain a positive density of su- 
perconducting electrons, suggesting that inhomogcncous 
LOFF phases are disfavored at small /. 

In some respects the same qualitative behavior we 
find here in the p-wave resembles what arose in s-wave 
[l^ . Namely, isotropic s-wave superconductivity has two 
branches of solution: the upper BCS which is stable and 
- for simple interactions - fully gapped, and the lower 
branch which has gapless modes but is unstable. The 
striking difference is that in p-wave the upper branch 
retains stability while developing a full two-dimensional 
fermi surface of gapless modes. Thus the anisotropic 
p-wave breached pair phase, with coexisting superfluid 
and normal components, is stable already for a wide 
range of parameters at weak coupling using the simplest 
(momentum-independent) interaction. This bodes well 
for its future experimental realization. 

In our model, which has no explicit spin degree of free- 
dom, gapless modes occur for either choice of order pa- 
rameter with residual continuous symmetry. By contrast, 
for '^He in the B phase the p-wave spin-triplet order pa- 
rameter is a 2 X 2 spin matrix, containing both polar 
and planar phases components, there are no zeros in the 
quasiparticle energies, and the phenomenolo gy broadly 
resembles that of a conventional s-wave state [l7| ! in the 
A phase (which arises only at T ^ ^|) the separate up 
and down spin components pair with themselves, in an 
orbital P-wave, and no possibility of a mismatch arises. 

Experimentally, the microscopic nature of the pair- 
ing state can be revealed most directly by probing the 
momentum distribution of the fermions, including angu- 
lar dependence. Time of flight images, obtained when 
trapped atoms are released from the trap and propagate 
freely, reflect this distribution. 

It is possible that the emergent fermi gas of gapless 
excitations develops, as a result of residual interactions, 
secondary condensations. Also, one may consider analo- 
gous possibilities for particle-hole, as opposed to particle- 
particle, pairing. In that context, deviations from nest- 
ing play the role that fermi surface mismatch plays in the 
particle-particle case. We are actively investigating these 
issues. 
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